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Our investigation is specially motivated by the stochastic version of a common 

model of potential spread in a dendritic tree. We do not assume the noise in the 

junction points to be Markovian. In fact, we allow for long-range dependence in 

Mh ' time of the stochastic perturbation. This leads to an abstract formulation in terms 

(-H I of a stochastic diffusion with dynamic boundary conditions, featuring fractional 

.4.^ , Brownian motion. We prove results on existence, uniqueness and asymptotics of 

Cy ■ weak and strong solutions to such a stochastic differential equation. 

Keywords: diffusion on network, fractional Brownian motion, strong solutions 
of infinite dimensional stochastic differential equations, invariant measures. 
MSC 2000: 34B45, 47D07, 60Hf5, 37L40 



o 
00 
o 

> 



1. Introduction 



> 

^— V \ This paper is concerned with a model of stochastically perturbed parabolic network 

'!;;j- ■ equations. We consider a network, i.e. a metric graph, on whose nodes we impose conti- 

nuity conditions, complemented by dynamic or Kirchhoff-type laws that also incorporate 
C3 ' stochastic noise. We employ known results concerning well-posedness and asymptotics of 

QP . the deterministic model associated with the equations in order to discuss solvability of the 

stochastic differential equation and convergence properties of the stochastic convolution 
process. 
^ , There is a well established theory concerning stochastic equations in infinite dimensional 

J^ \ spaces with additive Wiener noise, that can be found for instance in [11]. Our aim is 

to provide similar results for perturbation driven by a fractional Brownian motion. Since 
this noise is neither a semimartingale nor a Markov process, a different approach to a 
stochastic calculus with respect to it should be chosen. Our first result, Theorem 13.8^ 
provides the existence of a weak solution for our problem. In this case we extend the 
established results for Wiener noise. More interesting, and less classical, is the existence 
of strong solutions provided in Theorem 13.101 Even for the Wiener noise such a result, in 
the infinite dimensional case, is known only in some special cases; as a matter of fact, our 
assumption H > 3/4, together with the ambientation W{t) € D[{—A)°') for any a < 1/4, 
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which leads to if + a > 1, is the analog of the assumption W{t) G D{{—A)°^) for some a > 
1/2 required in [T] (since the Wiener case corresponds to H = 1/2). We shall also mention 
the paper [18], which inspired our proof, but where the stronger requirement W{t) £ D{A) 
is made. 

In Subsection 12. II we are going to introduce the deterministic model of the time evolution 
of potential in a dendritical tree according to Rail's linear cable theory. Successively, we 
are going to perturb it by stochastic terms that are represented by a fractional Brownian 
motion, as delineated in Subsection 12. 2[ Section [3] is devoted to the stochastic analysis for 
a network which contains only active nodes, while the general setting is sketched in Section 

m 

2. Biological motivations 

The motivation of this paper is provided by a stochastic model of diffusion for neuronal 
activity in passive dendritic fibers. The electrical behaviour of a neuron has been long 
studied, beginning with the pioneering experiments of H. von Helmholtz, more than 150 
years ago. In the 1950s, the ground-breaking papers of, among other, A.L. Hodgkin and 
A.F. Huxley [Hj and W. Rail [30] have represented a major breakthrough in the modern 
neurophysiology. Ever since, experimental results have dramatically increased. Successively, 
more and more interest has been revived by stochastic models - these may help in gaining 
a better understanding of neuronal behaviour and in predicting neuronal activity. 

A typical neuron cell consists of four parts: the dendrites, the synapses, the soma, and 
the axon. Dendrites are the input stage of a neuron and receive input from other neurons 
through axo- dendritic synaptical junctions. Synapses also connect dendrites to each other 
(leading to so-called dendro- dendritic interactions): in graph theoretical language, the den- 
dritic network is actually a tree converging into the soma, the cell's body. In [301 [3T] Rail 
introduced the so-called "lumped soma model" , representing (a part of) a dendritic tree as 
a lumped R-C circuit, thus allowing for an internal electrical activity of the soma. In Rail's 
model as well as in most of its generalizations (and in particular in leaky integrate- and- fire 
neuronal models, see [E]/ it is the soma's duty to elaborate these inputs and transmits 
them to other neurons through the axon. In this way Rail was able, under strong (geo- 
metric) symmetry and (electrotonic) isotropy assumptions, to describe the behaviour of a 
whole dendritic tree by a cable equation on a single "equivalent cylinder" (i.e., an individual 
interval of finite length) equipped with a dynamic boundary condition. Later on, multi- 
cylinder models have been proposed already in the 1970s in |32j-|33j in order to weaken the 
theoretical assumptions on the topology of the dendritic tree, thus describing more general 
classes of individual neurons. Such models have been extensively discussed in the 1990s, see 
e.g. [20] and subsequent papers. Aim of the present paper is to further generalize a model 
of Rail- type to a non-deterministic one, in order to allow for noise in the boundary, i.e., for 
non-deterministic inputs in the somata. 
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The starting point for the present analysis is represented by the recent article [26]. There, 
a thorough investigation of the deterministic version of such model has been performed: 
beside well-posedness results, several qualitative properties of solutions have been shown. 
Our mathematical approach also applies some of the methods of that paper. 

Rail's linear cable theory applies to the behaviour of the passive dendritic tree: the 
propagation of action potentials along an axon is modelled as a semilinear diffusion with 
a zero-order term that accounts for describing how the cell expends energy in order to 
propagate the signal, see below. On the other hand, it is known that in axons action 
potentials are initiated and propagated with fixed asymptotical speed and profile in the 
outward direction: an effect of regenerative self-excitation is thus needed while modelling 
axons. Such a phenomenon is commonly described by a coupled system of a semilinear 
diffusion equation and a nonlinear ordinary differential equation. Such a system has been 
introduced by R. FitzHugh and J. Nagumo as a simplification of the original Hodgkin- 
Huxley model: for a survey of this theory we refer, e.g., to [37]. Thus, the behaviour of 
the neuron in the dendritic tree and the axon is significantly different: as in [20] and [26] 
we will not attempt to introduce axons into the model considered in this note. A (naive) 
hybrid axo-dendritic model has been studied in a deterministic setting in [7]. 

Most of the aforementioned papers deal with a compartmental model of the neuron, 
subsequently investigated by numerical methods. Aim of this note is to provide a more 
theoretical approach based on an interplay of operator theory and stochastic analysis. 

Stochastic fluctuations of synaptic activity and post-synaptic elaboration of electronic 
potential is directly related to chaotic interferences in the streaming of charge carrying 
molecules. This chaotic streaming is classically modeled as a sequence of independent 
random variables; averaging on time, in dependence of the time scale, leads to choose a 
Brownian motion (or more generally a Levy process) as a model for the stochastic input 
process. Possibly combined with active properties of axons, this may lead to chaotic trigger- 
ing of action potentials. Several articles have actually dwelled on the stochastic dynamics of 
axons within the framework of the Hodgkin-Huxley or FitzHugh-Nagumo axonal models, 
see e.g. |38j or [36], mainly concerned with the issue of voltage fluctuations near threshold 
or the more theoretical approach in [4]. 

Less attention has seemingly been paid to the non-deterministic aspects of sub-threshold 
stochastic behaviour, either in passive or active fibers. However, there seem to be good 
reasons to perform such an analysis: in particular, "many computations putatively per- 
formed in the dendritic tree (coincidence detection, multiplication, synaptic integration and 
so on) occur in the sub-threshold regime" (quoted from [22] ) . A model for passive dendrites 
of infinite length with distributed as well as synaptic current noise has been proposed in 
jlT^llOj. whereas more experimental studies have been presented in [12]. A simplified model 
of a network with active nodes is proposed in [5]. Perhaps the clearest outcome of such 
investigations is that voltage noise increases with depolarization, cf. |16j . 
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Although our biological motivation leads us to consider a parabolic network equation, 
our results (with possibly the exception of the existence of strong solutions) are not confined 
to the 1-dimensional case. In fact, by means of the same techniques we may consider with 
minor changes also a stochastic initial-boundary value problem whose deterministic part 
is a diffusion equation with dynamic (a.k.a. Wentzell-Robin) boundary conditions on a 
bounded domain. The same remark applies to the stochastic analysis provided in Section 
[3] which could have been stated in abstract form and applied in different contests. 

2.1. Derivation of the equation. The basic object of our investigations, i.e., the dendritic 
network is modelled by a finite, connected graph G with m edges ei, . . . , e^ and n vertices 

Vi,...,V„. 

Let V denote the electric potential on the whole dendritic network, including the junctions 
(synapses and further ramification points) and the soma. Then V = V{t, x) is a function 
of the position x along the network and of time t only. Of course, V denotes the deviation 
from resting potential, which we rescale to (in concrete neurophysical measurements it is 
of approx. —IQmV). With a somewhat coarse but common approximation, we will assume 
the cell's body to be isopotential. 

By Vj we denote the electric potential on the edge e^, i.e., Vj{t,x) := V{t,Q.j{x)). Up 
to considering suitable rescaling diffusion parameters in the equations, we may and do 
parametrize the edges as intervals of unitary length. Actually, we avoid to introduce more 
general network elliptic operators: in fact, under uniform ellipticity assumptions this case is 
known to present no serious mathematical challenges over the basic case of a plain network 
Laplacian. 

Following Rail's passive cable theory, the transmission of post-synaptic potential in den- 
dritic trees can be mathematically described by a linear cable equation: on each edge ej 
(i.e., each minimal set of dendritic elements that can be reduced to an equivalent cylinder, 
cf. [20]), we consider the partial differential equation 

dV- d"^ 

(2.1) -^{t,x) = ^Vj{t,x)-pj{x)Vj{t,x), t>0, xG(0,l), 

We borrow from graph theory some basic notions, in order to describe the behaviour of 
the system at ramification points and to describe the architecture of the dendritic network 
in a compact form. We consider a graph G with vertex set V and edge set E through- 
out. We impose dynamical and time-independent conditions in the nodes vi, . . . ,v„(j and 
Vno+i) • • • ) v„, respectivelju. All nodes of the network are connected by the above intro- 
duced edges &j. We denote by J := J"*" — 3~ the nx m incidence matrix of the graph, where 

In biological applications a dendritic network is in fact a tree, in the graph theoretical sense, and therefore 
only one soma needs to be considered, i.e., no — 1. 
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3^ := (if,) and J~ := (i^) are the incoming and outgoing incidence matrices defined by 

■' 10 otherwise, ■' 10 otherwise, 

respectively. We also denote by r(vj) the set of all indices of those edges having an endpoint 
in Vj, i.e., 

r(vi) = {j G {1, . . . ,m} I 6^(0) = Vj or 8^(1) = Vj}. 

Motivated by Rail's lumped soma model (see e.g. |30j and [20j), we specify the behaviour 
of the system at the edges' boundaries by imposing several kinds of conditions in the nodes. 

• Both somata and synapses can be considered as isopotential. In other words, in all 
nodes Vj, i = 1, . . . , n, the potential V must satisfy the continuity assumption 

(2.2) Vj{t,Vi) = Vk{t,yi)=:qY for all j, fc G r(vi), t > 0. 

• In the somata vi, . . . , v^^, the potential qf {t) = V{t,yi) undergoes internal dynam- 
ics, subject to internal electrical activity and a stochastic feedback from the dendritic 
network. In the spirit of Rail's lumped soma model (see e.g. [20]) we impose gen- 
eral, possibly absorbing (and also possibly non-local) nodal conditions. They can 
be formulated as a Langevin-type equation 

(2.3) _gf(t) = _^i.^._^(t,v,)-5^6ih'?h'(i) + J^Qh^h(i), t>0, 

j=l h=l h=l 

where the last addend on the right hand side is a sum, over the edges incident in Vj, 
of (formal) derivatives of the stochastic inputs acting there. We shall give a precise 
definition and complete assumptions on the noisy terms Z in Section 12. 2[ 

• Motivated by the intrinsic randomness of external stimuli, we impose in the synapses 
and other ramification nodes Vj, z = tiq + 1, . . . ,ra, a version of a KirchhofF's law 
which is also perturbed by a stochastic term Zj. Such a condition can be equivalently 
formulated as 

m Qy n n 

(2.4) Y,tij-^{t,\/i) + Y,^ihqh =^CihZh{t), t>0. 

j=l h=l h=l 
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Summing up, the biological model we want to discuss is governed by the system of 
stochastic initial-boundary value problems 



(Vj{t,x) = V;'{t,x) - p,ix)V,it,x), 
Vj{t,Vi)=Ve{t,y^)=:qY, 



(2.5) 



t>o, xe (0,1), 

j = l,...,m, 

t>o, j,ier{yi), 

i = l,...,n. 



■y 



(t) 



J2 iijV'{t,yi) - Yl bihQh + S CihZh{t), t>0,i = l, 



,^0, 



j=l h=l h=l 

m n 



t > 0, i = no + 1, 



,n, 



X G (0,1), j = l,...,m, 



h=l j=l h=l 

V,{0,x) = Vo,{x), 
1/(0, vi) = qoi, 

For the sake of notational simplicity, let us introduce the Kirchhoff operators Ka,Kp 
mapping H^{0, 1;^") into C"° and C"'~'^°, respectively, defined by 



1, 



,no. 



I 



KaV:-- 






\ 



\ i=i / 



/ 



and KpV :- 



\ 



■ E ''no+ljK-(*) V„,j+i 
3=1 



■ E l^njV'{t,yr. 



A 



/ 



Thus, the vector KaV (resp., KpV) represents the differences between incoming and out- 
going flows in each of the active (resp., passive) nodes of the network. 
We introduce the matrices 



Ba 



"^h ;i=l,...,no' 



Br, 



hno+i)h}i=l,...,n-no ^nd B 



Br 



hi 



h I i=l,...,n 
h=l,...,n 



h=l,...,n h=l,...,n 

Roughly speaking, Ba and Bp encode the inhibitory and excitatory properties of active 
and passive nodes under the influence of the whole system. Also, in order to model the 
stochastic terms we introduce the matrix 



Ca 



h=l,...,no 



as well as 



i=l,...,no' ^P 
h=l,...,n 



'(no+i)/i jj=l,...,n-no 
h=l,...,n 



and C 



Ca 

Gri 



^ih I i=l,...,n 
h=l,...,n 



With no loss of generality, we may and do assume that the entries of the matrix C are non- 
negative real numbers and that there exist some entries in Ca which are strictly positive. 
In particular, the stochastic terms in ()2.3p - ()2.4p account for both the intrinsic cellular noise 
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and the external input signals. Observe that although the possibility of non-local conditions 
does not seem to be realistic in biological system, it may be interpreted as some form of 
external boundary feedback control of the system - e.g., with the aim of stabilization. 

2.2. Fractional Brownian motion. As customary in network applications, we use sto- 
chastic processes in order to model the input of our system. Altough the stochastic processes 
used in this field were often assumed to be Markovian, other considerations show that real 
inputs may exhibit long-range dependence, i.e., the behaviour of the process at time t does 
depend on the whole history up to time t. Another property that stochastic inputs show, 
also in telecommunication networks, is self- similarity: their behaviour is stochastically the 
same, up to a space scaling, on changing the time scale. These features are common in 
physiological systems and have been already applied to the neuroscience, see e.g. [27], [9], 
and the contributions in [341 Part II]. Among the processes which exhibit these proper- 
ties, we propose to model our stochastic input with fractional Brownian motion (fBm for 
short). This class of processes also verify stationarity of the increments and continuity of 
trajectories. 

A 1-dimensional fractional Brownian motion is a centered Gaussian process {B^{t), t > 
0}, such that 

E[S^(t)S^(5)] = -^(i'"" + s^"" -\t- sD. 

2x3 

The constant H is called Hurst parameter and takes value in (0, 1). The construction of 
fBm was proposed by Mandelbrot and van Ness [21] using the representation 

JQ J-oo 

where B = {Bt, t > 0} is a standard Brownian motion. We shall be interested in a different 
characterization of fBm, based on the observation that, given the fBm B (t), there exists 
a unique Brownian motion W{t), adapted to the same filtration, and a kernel Knit, s) such 
that the identity 

5^(t)= / KH{t,s)dW{s). 
Jo 

holds. In case H = i^, B^'^{t) = W{t) is a standard Brownian motion. In this paper we 

restrict our interest to a particular class of fBm and impose the following. 

Assumption 2.1. We assume that there exists a probability space (fi, 9", P) endowed with a 
filtration {3't, t > 0} that satisfies the standard assumptions, i.e., it contains all negligible 
subsets of i7 and it is right- continuous. 

On this space we are given a family of independent standard Brownian motions {Wh{t), t > 
0}/i=i,...,n; then we consider, for each h = 1, . . . ,n, the fBm 



Zh{t)= I KH{t,s)dWh{s), t>0. 

JO 
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We let Z"'{t) = (Zi(t), . . . , Z„g(t)) and Z'^{t) = (Z„,j+i(t), . . . , Z„(t)) denote the stochas- 
tic inputs in active and passive nodes, respectively, and Z[t) = {Z^{t),ZP{t)) be the n- 
dimensional stochastic input process; we refer to Z{t) as a n- dimensional fBm. 
The Hurst parameter H belongs to (2; !)• 

2.3. Wiener integrals for fBm. For any given Hilbert space H, we are interested in 
the introduction of Wiener-type integrals for H-valued, square integrable functions, defined 
with respect to a fractional Brownian motion Z{t). This is a known topic in literature and 
we follow the stochastic calculus of variations approach, see the monograph [2] and the 
references therein; we shall sketch below the notation and the main results. 

In the Brownian case ff = i, the space of Gaussian random variables defined as .^ = 
J^f{s)(iWs is isomorphic to the space L^(]R;H). This isometry is classically proved first 
for step functions and then is extended to the whole space. A similar construction was 
first investigated in [29] for the fractional Brownian motion. Let £(H) be the set of step 
functions $ : R ^- H 

m 
i=0 

the stochastic integral /(^) is defined by setting 

/($)= / ^s)dB^{s) = Y,MB^{t,+^) - B^iti)) 

and it is a Gaussian random variable with zero mean and covariance that can be computed 
using the covariance matrix of B . The set of Gaussian random variables defined by the 
elements of £(S) is a subset of the first Wiener chaos 

:Ki = {XGL\n-E)\3{fn)n>oC8.{E) : lim /(/„)= X in L^{n-E)} 

n— >oo 

generated by B . The reproducing kernel Hilbert space (RKHS for short) A is defined as 
the closure of £(H) with respect to the scalar product 

(l[o,t]|l[o,s]) = R{t,s); 

the mapping / : Ifq^^i — > B^{t) defines an isometry between A and "Ki. 

As opposite to the Brownian motion case, for H > -^ the RKHS A for the fBm B (t) 
cannot be identified with a space of functions; however, A contains (linear subspaces which 
are isometric to) inner product spaces of functions. Proceeding as in [29] we can prove that 
this space can be identified with the space |!K| of functions in L^(]R; H) n L^(]R; H) endowed 
with the norm 

U=»H [ [ |t-rp^-2|</>(r)|s|</'(t)|Hdrdi. 
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This space contains the class of step functions £(H) and the isonietry 



E 



(/.(t)dS^(t) 



2 

1^1 



holds. 



3. The abstract formulation 

Before discussing the complete stochastic differential equation in Section HI we begin by 
considering the system with the following simplified form of (12.4^ : 



(3.1) 



m rjT^ n 

j=l h=l 



hql = 0, 



t >o. 



To begin with, we introduce the Hilbert spaces X := (L^(0, 1))™ and dXa := C"". On the 
domain 



D{A) :-- 



d := 



V 



^V I G (^'(0, 1))"' X C"o s. th. 3g^ G C" with (J+) ' g^ = V{0), 

(r)Tg^ = y(l), {qY,...,ql) = {q^„...,qLo)^ and KpV = Epq"" 
we define the operator A by 



(3.2) 



Ao:=A{ V 



V" - pV 
KaV - Baq^ ' ' 



where V" denotes the vector (T^/', . . . , T^) G X consisting of (one-dimensional) second 
derivatives of the real- valued entries of the vector- valued function V = (Vi,...,Fm) G 
(/f'2(o,l))'". 

One can directly check that the initial value problem associated with (|2.ip - (j2.2p ~ (|2.3p - 
(|3.ip can be equivalently formulated as an abstract stochastic Cauchy problem 



(3.3) 



do(t) = AX){t) dt + Qa d^"(i), 
0(0) = Do, 



t> 0, 



where the initial value is given by do := (Vo, qo)^ 6 X := X x dXa and Qa maps every vector 
q £ dXa in (0, Caaq)^ G X. 

Our first aim is to collect some results on the underlying deterministic model, which will 
show well-posedness and further qualitative properties of our system. Although formally 
new, the results of this section can essentially be proved combining the ideas of |24j and [26], 
where parabolic network equations with passive only, non-local interactions and local, active 
node conditions have been considered, respectively. 
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In order to prove the generation property of the operator A, we mtroduce a Hilbert space 



V 



\ G (i/^O, l))*" X C"» s. th. 3 r?^ G C" with O+Yq^ = T/(0), 



(r)^g^ = v(i), (^r, ■ ■ ■ ,'?ro) = (9^1, • • • ,C) 



and the sesquihnear form o : V x V ^ C defined by 

a(u, tt) := ao(u, o) + ai(u, o) := {U'\V')x + {{pU\V)x + {Bq^\q^)dx^) ■ 

We emphasize that this form is in general not symmetric, but it is always densely defined. 
Mimicking the proofs of |241 Lemma 3.4] and |26l Lemma 3.3] one can prove that the 
operator associated with a is A. Observe that o and hence A are self-adjoint if and only 
if the scalar matrix B is hermitian, i.e., if and only if mutual interaction of the nodes are 
symmetric. 

Combining the techniques presented in [26] and [2l] for discussing dynamic and non-local 
boundary conditions, respectively, we can prove the following. 

Proposition 3.1. Let B G Af„(C) and p G (L-'^(0, 1))™. Then the operator A generates 
a strongly continuous, analytic and compact semigroup (S(t))i>o on the Hilbert space X : = 
X X dXa. 

Proof. We first observe that the leading term in the form a, i.e., Oq, is clearly X-elliptic and 
continuous. Furthermore, it follows from the Gagliardo-Nirenberg inequality that the space 



C(G) : = 



\ ) G (C([0, 1]))™ X C"" s. th. 3 g^ G C" with {'i+Yq^ = F(0), 



(r) V = v{i), (^r, ...,ql,) = {qL ■ ■ -^qLo) 



V G (C([o, 1]))"* s. th. 3q^ eC" with (a+)^g^ = v{o), (r )"^Q^ = y(i)} 

of continuous functions over the network is embedded in an interpolation space of order ^ 
between V and X, and obviously ai : C(G) x C(G) — *■ C is bounded. Thus, by a suitable 
perturbation argument, cf. [25l Lemma 2.1], one concludes that also their sum a is X-elliptic 
and continuous. Thus, by [281 Prop. 1.51, Thm. 1.52] the associated operator A generates a 
strongly continuous, analytic semigroup. Moreover, by Rellich-Khondrakov's theorem the 
imbedding V ^-> X is compact and accordingly the semigroup is compact. D 



Remark 3.2. If coefficients bij and Pj{x) are real constant, resp. real valued functions, 
then S(t) maps real valued functions in real valued functions, thus, despite of the choice 
of complex valued function spaces, we may appeal to this remark in order to justify the 
application of our results to biological models. 
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All further properties of the semigroup essentially depend on the matrix B: e.g., it is 
possible to characterize asymptotics of the semigroup generated by A by means of positivity 
of the L^-coefficient p and/or of the definiteness of the (in general non-hermitian) matrix 
S, cf. [El Thm. 2.3]. 

Proposition 3.3. If p > and B is positive semidefinite, then the following assertions 
hold. 

(1) (S(t))t>o is contractive, i.e., ||S(i)|| < 1 for all t > 0, and hence mean ergodic. 

(2) If additionally p = and B*l ^ 0, or else if B = and p ^ 0, then (S(t))t>o is 
strongly stable, i.e., limt^oo §(i)ti = for all D G X. 

(3) If additionally p > po > and/or B is positive definite, then (S{t))t>o is uniformly 
exponentially stable, i.e. ||S(t)|| < e""^* for some a; > and all t > 0. 

By compactness of (S(i))j>o and an abstract criterion due to Ouhabaz, see [281 Thm. 
2.6] we obtain the following. We do not sketch here the proof and refer the reader to [26] 
and [24| for technical details. 



Corollary 3.4. Ifp = and B = 0, then {§{t))t>o converges towards the (strictly positive) 
projection onto the 1- dimensional eigenspace spanned by the constant vector 1 £ X x dXa. 

Assume for a moment that A is negative definite. Then, since yi is a sectorial operator, 
we can introduce the scale of interpolation spaces defined with respect to the domain of 
the operator A. The leading term of the form associated with A is symmetric, i.e., A is 
self-adjoint up to dropping the boundary conditions given by B, which are only relevant for 
higher powers of A. 

More precisely, it follows by the spectral theorem that complex interpolation spaces 
Xq, := [D{—A),X\a are equivalent to the domains of fractionary powers {—A)'^ of —A. 

Factoring A as in [231 Lemma 2.7] one can (up to similarity) decouple its domain into 
the product IIq{G) x dXa. Here IIq{G) denotes the space of {II^{0, l))"*-functions with 
Dirichlet boundary conditions in the active nodes vi , . . . , v„,(, . For a < ^ , it follows that 
D{{-A°') is isomorphic to ij2"(G) x dXa. 

In general, it is always possible to perform the same steps by defining the complex 
interpolation spaces Xq, := [D{—XI — A), X]a for a suitable constant A. 

3.1. Stochastic differential equations with additive fBm. In this section we solve the 
stochastic differential equation (|3.3p 



fdt)(i) =yid(t)dt + eadZ"(t), t>0, 
''■'* 1 .(0) = .„, 

where Oq := (Vq, qo)'^ £ X := X x dXa. 
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Definition 3.5. Following the Wiener case treated for instance in [DPZ92] we call an X- 
valued adapted process d = {D(i), t € [0,T]} a strong solution to \3/J\) if V has a version 
such that 

(51) for almost all t G [0,r], P(d(t) G D{A)) = 1; 

(52) for any t G [0,r], P ( / \A\){s)\\ds < +oo J =1, and 

(53) for any t G [0, T] there holds 

(3.5) 0(t) = Oo+ f Av{s) ds + Ca Z^it). 

Jo 

Definition 3.6. Similarly, we define a weak solution to (3.3\) if 

(Wl) for any t G [0,T], P ( / |d(s)||ds < +oo J = 1, and 
(W2) for any t G [0,T] and any y G D{A*) it holds 

(3.6) (d(t),y) = (do,y) + (/ v{s)ds,A*y) + {eaZ^{t),y). 

Jo 

3.2. Existence of weak solutions. In order to provide the existence of a weak solution, 
our main interest lies in the stochastic convolution process 

WA{t):= f S(i-a)e,dZ». 
Jo 

The key point is to find estimates for the variance, like 

no 
nWA{t)\l2 = 5^ |||§(t - a)eaek\^, l(o,t)(cT) 



|2 

Ia 
fc=l 



where A is the RKHS associated with the fBm Z". This means that in case H = 1/2 (the 
Wiener case) A is the space L'^(R) and in case H > 1/2 we can take instead of A the space 

1^1- 

Since §(t) is a semigroup, there exist constants M > 1 and w G M such that ||§(i)||L(X) ^ 
Me"^*; however. Proposition 13.11 implies that we can choose M = 1 and u> < 0. It follows 
that 

nwAit)\l2 < M\\ea\\Hs ||e"" i(o,t)(^)|ll ; 

in the Wiener case, the norm is given by L e^^^ ds = 2Zj(^'^* — 1) < W~V °^ ^^^ other hand. 



in the fBm case, the estimate becomes 

(3.7) / I e'''e''''\s-rf^-'^dsds<Const.-^\\-e 

Jo Jo ^^^ 
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and we get 

(3.8) sup E|MOi(i)||< Const. \\Qa\\]is- 

iG[0,T] 

With a little effort, it is possible to be more precise on the regularity (in space) of 
the stochastic convolution process. We remark that the case H = 1/2 is treated in \\\\ 
Proposition 6.17]. 

Lemma 3.7. Assume that the operator A is injective and dissipative. Hence, we can 
consider the fractional powers of the operator A and the family of spaces Xq, with equivalent 
norm \\v\\x^ := \{-A°'v\x. 

Let H > 1/2 and assume that 

H-a>0. 
Then the stochastic convolution process verifies 

(3.9) E f \\WAt)fx^dt< Const. 

Jo 

Proof. We get 

E/ \\WAit)\l^dt<Y, ||||S(t-a)eaefc||^^l(o,i)(a)||^dt 
Jo i^^^ Jo 

By definition of the norm in Xq, we may compute 

||||S(t-a)eaefc||x„ l(o,i)(a)||' < [ [ \{-ArSit-a)eaek\\i-ArS{t-e)eaek\\c7-e\^''~^dade 

Jo Jo 

-^ ^ Jo Jo ^"^" 
<M^Ca,H ICaefcl^t^-" 

and for any finite T > the quantity in (j3.9p is finite. D 



Theorem 3.8. If A generates a Co-semigroup (§(i))t>Oi then there exists a unique weak 
solution to 113. 3\) and it has the representation 

(3.10) t)(t) = S(t)t)o+ / S(t-s)eadZ"(s), te[0,T]. 

Jo 

Proof. Fix t G [0,T]. Let us first notice the following identity which holds for every ^ G 
C\[0,t];D{A*)): 

(3.11) (t)(t),cl>(t))-(do,$(0)) 

ft rt ft 



(d(s),6(s))ds+ f {v{s),A*^s))ds+ I («5(s),eadZ"(s)). 

Jo Jo 
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Assume now that u is a weak solution to ()3.3p . Choosing <I>(s) = S*(t — s)yl(o,t)('S) in ()3.1ip 
we get 

(D(t),y)-(Do,§*(t)y) 

= / (0(5),^S*(t-s)y)ds+ / {x,{s),A*2>*{t-s)y)ds+ [ (S*(t - s)y, C, dZ«(s)) 
Jo Jo Jo 

that is 

(d(t), y) - (S(t)t)o, y) = / (y, S(t - s)e, dZ«(s)) 

JO 

which imphes uniqueness of the weak solution. 

We now proceed to show that the process u defined in (IS.lOp is a weak solution to p.3p . 
Taking in mind formula (j3.6p we compute 

/(§(cj)Do+ / §ia-s)edZ''{s),A*y)da 
Jo Jo 

= / (do,S(a)yi*y)dcT+ / [ {S{a-s)eadZ^{s),A*y)da 

Jo Jo Js 

= f^{^,,§*(a)y)da+ f\ f ' As*(a)y da , C, dZ»(.)) 
Jo Jo Jo 

= (S(t)Oo, y) - (Oo, y)+ f {§>*{t - s)y - y, C, dZ''{s)) 

Jo 

= (S(i)oo, y) - (do, y) + / {y, S(t - s)e dZ"(s)) - e, Z'^(t) 

Jo 

and the representation (|3.1Up yields to the result. D 

Corollary 3.9. Assume that A is a bounded linear operator. Then the stochastic convo- 
lution process 

(3.12) WA(.t):= I S(t-s)e„dZ'^(s) 

Jo 

verifies 

(3.13) WA{t) = [ AWa{s) ds + Ca Z''{t). 

Jo 

3.3. Existence of a strong solution. 

Theorem 3.10. Let A be a sectorial operator and assume that Ca rnaps X into Xq for 
any a < 1/4. Assume that the Hurst parameter H of the fBm Z°-(t) verifies H > 3/4 and 
take the initial condition do € D(A). Then there exists a unique strong solution to equation 
[KM). 
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Up to rescaling, since we deal with a finite time interval [0,T], we can assume that A is 
dissipative and injective. 

Yosida approximations. Before to proceed to prove the existence of a strong solution for 
()3.3p . we need to introduce the family of operators yi„, n G N, defined by 

An = nAR{n,A) = n^R{n,A) - nl. 

This family has been introduced in Yosida's own proof of the celebrated Hille- Yosida gen- 
eration theorem. It enjoys remarkable properties: the operators An are bounded; they 
commute with one another as well as with A and with the resolvent operators of A; and 
finally 

lim 2>n{t)x = ${t)x for all x S X and lim AnX = Ax for all x E D{A), 

n— >oo n— ►oo 

by [131 LemmaII.3.4]. 

For future references we need a lemma concerning the speed of convergence of the semi- 
groups generated from Yosida approximations to S(t). 

Lemma 3.11. Let A he a sectorial operator. Let also C G L(X, Xq) for some a G (0,1). 
Define, for any n G N, the sequence of functions fn{(^) setting 

/„(a) = |K(S„(a)-S(<T))e||i(X). 
Then for given T > it holds 

lim sup |/„(a)| = 0. 
"^°°(7e(o,T) 

Proof. Write 

\\Ani§nia) - S(a))e||i(x) = n\\A^-^Rin,A)i§nia) - S(a))yi"e||i(x) 

<n\\A'-^Rin,A)\\L(x)\\iM<^)-H<T))hm ll^"e||i(x) 

We estimate, for cv = s{A) the spectral bound of A and some e > 0, for any n > to + e, 

r-oo 

\\A^-^R{n,A)\\L(r)< / e-"*||S(tMi-"||z.(x)dt 
Jo 

/>oo 

Jo 

We now recall that in [3] it is proved that \\R{\An) - R{\A)\\ < ^i±^. Define S(t) by 
the Dunford- Taylor integral 



§(a) = — [ e^''R{X,A)dX. 
2-1X1 Jy 
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Then we get 



{%^{a)-%{a))\\Lix) 



< 



-^ [ e^''[R{X,A)-R{\,An)]dX 



L{X) 



where F is a piecewise smooth curve consisting of three pieces: a segment Fi = {re"*^^"'^' : 
1 < r < oo}, a segment F3 = \re^^^~''> : 1 < r < cxd} and the arc F2 = {e*^ : — (^ — e) < 
/? < (^ — e)}. Take the norm inside the integral; it follows 

1 



(^) - S(^))IIl(X) < ^ y^ le""^! I|[^(A,^) - R{X,An)\\\L(^) dA 



<^i±^/|e-|dA 
invr 



(1 + Mi)2 
2n7r 



|,-|dA<(l±^ 



n 



U 



Proof of Theorem 13.101 It is possible to take the initial condition Oq = 0, possibly 
considering the process tt(t) = o(t) — §(t)Do- Then we shall prove that t)(t) = Wji{t) is a 
strong solution of (j3.3p with zero initial condition. 
We notice that, by Corollary 13.91 the identity 



I S„(t-s)e,dZ'^(s)= I AnWA„{s)ds + QaZ' 

Jo Jo 



Wju.it) 



holds. We shall now send n to infinity in order to get ()3.5p . 
As a first step we prove that 



(t). 



(3.14) 
Notice that 



lim sup E\WAAt) - W^im = 0. 



WA^t) - WA{t) = [ [Sn(t -S)- §{t - s)]eadZ''{s) 

Jo 



and setting q/. := Ce^ (where {e^} denotes the canonical basis of dXa) one concludes that 

Ix 



nWA^it)-WAit)\%- = E 



V / [§n{t-s)-§it-s)]qkdZkis] 



X 



no 



J]|||[S„(.)-S(-)telxl 



(o,t)l-JII|:K| 



fc=i 



Hence, since it holds that ip{s) := [Sn(s) — S(s)]efc|x converges to as n — > 00, uniformly in 
s e [0, T] and 



l^(-)l{o,t)(-)ll|:K| 



t Pt rp2H 

/ \i;{r)ms)\\s-r\^''~'dsdr< sup 

Jo H[ZH - ij selO,T] 
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the claim ([31^ holds. 

Next we show that Wyi(t) G D{A), P-a.s., for all t e [0,r]. In the case of a Q- Wiener 
process W{t), the thesis follows, for instance, from an application of [DPZ92: Proposition 
4.15]. In our case, the assumptions can be modified: in fact, we require 

no 

(3.15) Y. ||l'^S(^)eaefe|xl(o,t)('^)|||:K| < +°°- 

fc=i 

We can then mimic the proof of the aforementioned result by Da Prato-Zabczyk and obtain 
the claimed result. Just observe that Qa^k belongs to Xq, for any a < |, hence 

|yiS(o-)eaefc|x < Const. r-\ t € [0,r], a < 1/4. 

A direct computation shows that 

(3.16) r Tr^-^s-V- sP^'-'drds = Const. T^^+^''-^ 
Jo Jo 

provided a -\- H — 1 > 0; notice that the existence of an a which satisfies the bounds 

1 - iJ < Q < 1/4 only holds if i^ > 3/4. 

Now, our last claim is to prove the following convergence 

(3.17) lim sup E\AnWA„{t) " •AWOi(t)lx = 0. 

Together with (j3.14p . this will imply that we can pass to the limit in (j3.13p and prove p.5p . 
hence the claim follows. 

We split the estimate in two parts, since 

(3.18) ElAnWA^t) - ^Wjiml < 2(^n,l(i) + ^n,2(t)) 

with 

^n,l{t) = ElAnWA^t) - ^nWA{t)\l 
$„,2(t) = nAnWAit) - AWA{t)\l. 

Notice that 

E\AnWAAt) - -^nWA{t)\l < C, \\\An (S„(t - s) - §{t - s)) e«|i.(x) || |^| 

and this term was treated in Lemma |3. lit it holds that sup^gj-Q^^ /«(s) ^ as ?i — > oo and 
from Lebesgue dominated convergence theorem it follows that ^n,i{t) -^ 0. 

Next we proceed with the second term in ()3.18p . We can follow the same steps to get 

no 

II I r /T A ic/j. „\ ro „_ \ _ J, 

\\Ji\ 



$n,2(t) = E\AnWAit) - AWAit)\l = Y\\\[A- An]§it - s) e,efc|x|'' 



fc=l 
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Now since [yi— yi„]x ^ as n — > oo for any x € D{A), it is sufficient to show that each term 
is uniformly bounded by an integrable function and apply dominated convergence theorem 
in order to obtain the desired convergence. Further, recall that An is a bounded operator; 
hence, proceeding as above, we get the bound (compare with (13.160 ) 

II II 2 

|||[yi -yin]S(s) eaefc|xl(0,t)(s)|||^| < +00. 

In conclusion we obtain that both <I>„,^i(s) and $n,2(s) converge to as n — > oo uniformly 
in s G [0, T] and the last claim is proved. 

3.4. Long time behaviour. In this subsection we will be concerned with the asymptotic 
behavior of Wji{t) for t — > oo. Let us recall from |1H Theorem 11.7] the equivalence between 
the existence of an invariant measure and the boundedness of the covariance operator Qt 
for Wj{{t) at all times i > 0. In the truly dissipative case we obtain the following. 

Proposition 3.12. Assume that p > po > and/or B is positive definite. Then Wji{t) 
converges as t ^> oo to a centered Gaussian random variable p.. Furthermore, fi is an 
invariant measure for equation 113. 3\) . 



Proof. The assertion follows if we prove that 

sup TrQt < oo, 
t>o 

where Qt is the covariance operator of Wji{t). We have discussed this quantity in previous 
section and we know that 

no 

nWA{t)\Jc = E |||S(^)eaefc|xl(o,t)(^)||l 
fc=l 

Now since in our assumptions, by Proposition 13.31 (3) the semigroup is uniformly exponen- 
tially stable, there exist M > 1 and lo > such that 



/•oo 

Tr Qt<M\\ea\\HS / e-2'^'" da < Const. 
Jo 



In case Zh is a fBm, we use in the space A the norm \'K\ and we get (compare (|3.7p ). for t 
large. 



no 



nWA{t)\l < M^ |||S(a)eaefc|x||f^| < Const. \\Ca\\HS 

k=l 

This concludes the proof. D 



1 t2H-2 
+ T- 



CJ2^ 



+0 as t-^oo 
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We study now the uniqueness of the invariant measure and we make precise the weak 
convergence of the associated transition probabihties. Notice that we can appeal again to 
the theoretical results in [HI Ch. 11]; hence, our results are based on the properties of the 
semigroup S(t) obtained before. 

Proposition 3.13. Assume that p > and B is positive semidefinite; then the following 
assertions hold. 

(1) If additionally p = and B*l ^ 0, or else if B = and p ^ 0, then there exists at 
most one invariant measure fi for i3.3\) . 

(2) If additionally p > po > and/or B is positive definite, then there exists one 
invariant measure fj, for Ii3.3\) . Further, for all initial condition Oq the solution d(t) 
converges in law and in L^(0; X^) to fi. 

In different terms, the situation that we obtain is that there exists an equilibrium state 
(which can be thought of as resulting in the long term behavior) for the neuronal network. It 
may be interesting, in this connection, to study large deviations to estimate the probability 
of onsets of chaotic impulses, compare [35|[5]. 

The situation is different in case there is no dissipation acting in the network. Then by 
Corollary 13.41 the semigroup converges toward a projection and for some to there holds 



\\§maek\\>hl\xY,\{lXaek)\ t>to 



k 

which is bounded below by the sum of the entries of k-th. column of the matrix Ca- Since 
at least one entry of Ca is strictly positive, as t goes to infinity it follows 

no 

^|||S(a)eaefc|l(o,i)((T)|||^l > Const. ||l(,„,i)(a)|||^| > Const, (t - tof" . 

k=l 

Therefore, the trace of the covariance operator Qt becomes unbounded as t ^ oo indipen- 
dently of the kind of noise that we consider; this (negative) result is interesting enough to 
be recorded in a proposition. 

Corollary 3.14. In case no dissipation occurs in the network (p = and B = 0) then there 
does not exist an invariant probability measure for the system. 

4. General stochastic perturbation 

In this section we briefly discuss a possible approach to inhomogeneous boundary value 
problem. This method has been proposed in [19] and has been extended to stochastic 
differential equations in [6]. 
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Let dX"^ := C"~"° and define a boundary operator R : D{A) — > dXp by 

R^ = R{V,q^f:=KpV-Bpq^. 
Consider the space 

V := ( ^y] G D{A) X dX^ s. tli. 3 g^ G C" with {3+)^q^ = 1/(0), 



D{A) :-- 



(r)V = V^(l), {qY,--.,ql) = {qL---,Qano), and i?0 = 0^ 



inoy Wal) ■ ■ ■ 1 Hano 

and the operator on X"^ x 5X^ x dX^ defined by 

'A\} 



Av.-, Q 

Then, it is known (cf. [191 § 3]) that A generates an analytic semigroup because A does 
so. Such a semigroup is never stable, since it acts as the identity on the second component 
of the vectors in (X^ x dX^) x dXp. In fact, even if (S(t))i>o is uniformly exponentially 
stable, the semigroup generated by A is given by 

g(i)^= /"§(*) (/-§(t))D^''' 

where Dq ' is a so-called abstract Dirichlet operator, cf. [191 Lemma 3.3]. 

We can now examine the case of a system featuring the presence of (stochastic) inputs 
in the passive nodes. The system is described by the operator matrix A on the space 
(X^ x dX^) X dXp. Let Z{t) be the n-dimensional stochastic process which models the 
input in the (active and passive) nodes, and C = (0 Ca Cp)'^ be the covariance operator of 
Z{t), where Ca and Cp are the matrices defined in Section 2. Then the stochastic model 
can be written in the form 

du{t) =Au{t)dt + CdZ{t) 
^^■^^ u(0) = no 

which is solved in mild form by 

(4.2) u{t)=S{t)uo+ / S{t-s)CdZ{s). 

Jo 

Proceeding as in Section [3. II we obtain the following result. 

Theorem 4.1. Assume that uq G {X'^xdX'^)xdXp. Then the process u = {u{t), t G [0, T]} 
is a weak solution of i4-l\ )- 

Assume further that uo G -D(A) and H > 3/4. Then the weak solution is a strong solution 
ofgJl). 
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We can be more precise in formula ()4.2p since we know the representation of the semigroup 
S(i): 

'S(t) (/ - S(i))D^'^' 
hence 



^(^) Mo 



u{t) 



Ua{t) 
CpZ{t) 



where 



Ua{t) = [ Ht- s) Qa dZ{s) +/"(/- S(i - s))V^'^Cp dZ{s). 
Jo Jo 

In this case, whenever the matrix Cp is not identically zero, we are again in presence of a 

gaussian process with finite trace class covariance operator at any time but not bounded as 

i — 5- oo and, similarly to Corollary 13.141 we obtain the following result. 

Corollary 4.2. If the behaviour of passive nodes is affected by some stochastic inputs (i.e., 
the matrix Cp is not identically zero), then there does not exist an invariant probability 
measure for the system. 
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